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Abstract. 

We  show  how  the  study  of  variations  in  orientation  of  a  terrestrial 
reference  system  (TRS)  in  space  may  be  done  directly  in  terms  of  the 
motion  of  the  pole  of  the  TRS  and  rotation  around  it,  and  how  a  separa¬ 
tion  of  these  variations  into  low  frequency  and  high  frequency  (retrograde 
and  prograde  diurnal,  semidiurnal,  •  •  •)  bands  enables  one  to  character¬ 
ize  and  model  variations  belonging  to  the  various  bands  and  to  estimate 
them  simultaneously  from  observational  data  by  a  uniform  procedure. 
Introduction  of  the  Celestial  Ephemeris  Pole  (CEP)  or  other  Celestial 
Intermediate  Pole  (IP)  is  not  only  unnecessary,  but  also  gives  rise  to 
needless  debate  as  to  whether  variations  due  to  particular  causes  are  to 
be  included  in  the  celestial  motion  of  the  IP  or  in  its  terrestrial  mo¬ 
tion,  and  leaves  the  question  of  estimation  of  high-frequency  signals  in 
either  frame  unresolved.  In  regard  to  UT1,  we  point  out  that  the  “cor¬ 
rection  terms”  through  which  the  concept  of  the  nonrotating  origin  is 
implemented  emerge  naturally  from  fundamental  kinematical  relations, 
and  use  this  observation  to  identify  the  correction  terms  to  be  employed 
when  the  Earth  orientation  parameters  are  defined  in  relation  to  the  pole 
of  the  TRS  rather  than  an  IP. 


1.  Introduction 

It  has  been  a  widely  held  belief  that  estimation  of  the  high-frequency  part  (with 
periods  below  2  days)  of  the  variations  in  Earth  orientation  cannot  be  done 
by  analysis  of  VLBI  data  gathered  routinely  in  one-day  sessions  at  intervals  of 
several  (typically  5  to  7)  days.  The  notion  of  a  Celestial  Ephemeris  Pole  (CEP) 
which  has  no  high-frequency  components  to  its  motion  in  terrestrial  as  well  as 
celestial  reference  frames  is  a  counterpart  of  this  belief  on  the  conceptual  front, 
though,  a  causal  connection  between  the  two  is  hard  to  pinpoint.  Both  of  these 
ideas  are  reflected  in  the  traditional  (and  current)  procedure  for  analysis  of  VLBI 
data,  which  decomposes  Earth  orientation  variations  into  separate  motions  of 
the  CEP  in  the  celestial  and  terrestrial  reference  frames,  with  no  high-frequency 
content  in  either,  and  a  nearly  diurnal  rotation  around  this  pole.  The  investment 
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made  over  a  couple  of  decades  in  algorithms  and  software  which  implement  this 
procedure  is  one  of  the  potent  factors  which  deters  any  significant  move  away 
from  this  pole,  and  understandably  so. 

Yet,  higher  frequencies  are  present  in  Earth  rotation  variations,  and  their 
amplitudes  can,  in  fact,  be  estimated  with  considerable  precision  from  existing 
VLBI  data  sets,  notwithstanding  the  belief  mentioned  at  the  outset:  Herring 
and  Dong  (1994)  showed  how  when  they  estimated  the  amplitudes  of  ocean-tide- 
induced  variations  having  diurnal  and  semidurnal  frequencies  in  the  terrestrial 
frame.  Recent  works  on  rigid  Earth  nutations  have  shown  that  diurnal  and 
semidiurnal  frequencies  occur  also  in  the  motions  of  the  Earth’s  axis  in  space. 
The  concept  of  a  slow-moving  CEP  has  thus  lost  its  relevance. 

The  main  question  addressed  here  is:  do  we  need  a  replacement  for  the 
CEP?  To  put  it  differently,  would  an  Intermediate  Reference  Pole  help  us  to 
accomplish  something  that  cannot  be  accomplished  without  it?  We  believe  the 
answer  to  be  in  the  negative:  an  intermediate  pole  is  redundant.  We  outline  our 
reasoning  here. 

Stated  briefly:  The  information  on  Earth  orientation  provided  by  ground- 
based  observations  of  celestial  objects  relates  to  the  orientation  of  an  Earth- 
fixed  reference  frame,  and  in  particular,  that  of  its  third  axis.  Even  when  Earth 
orientation  parameters  (EOP)  are  defined  in  relation  to  the  CEP,  their  variations 
have  to  be  inferred  from  the  observable  quantities,  which  relate  only  to  the 
motion  of  the  TRS  in  space.  Even  with  the  present  methods  of  VLBI  data 
analysis,  they  are,  in  actual  fact,  obtained  in  this  manner,  as  will  be  made 
manifest  below.  The  seeming  role  of  the  CEP  merely  masks  this  basic  fact. 

The  definition  of  the  axial  rotation  of  the  Earth,  and  hence  that  of  UT1, 
is  tied  to  the  choice  of  the  axis  about  which  the  rotation  is  to  be  defined.  The 
concept  of  the  nonrotating  origin  (Guinot,  1979)  may  be  implemented  with  any 
chosen  axis.  We  draw  attention  here  to  the  fact  that  the  small  “corrections” 
(Capitaine  et  al. ,  1986)  to  be  applied  to  obtain  UT1  from  the  finite  angle  of 
rotation  about  that  axis  —  defined  in  the  decomposition  of  Earth  rotation  into 
motions  of  the  pole  and  axial  rotation  around  the  pole  —  follow  directly  from 
the  kinematical  relations  connecting  the  Earth  rotation  parameters  (spin  rate 
variation,  and  wobble)  to  the  Earth  orientation  parameters. 


2.  Transformation  from  CRS  to  TRS  via  an  intermediate  pole 

The  transformation  from  a  celestial  reference  system  (CRS)  to  a  terrestrial  one 
(TRS)  via  an  intermediate  pole  (IP)  is 

T  =  WRS,  (la) 

where  S  takes  the  pole  of  the  CRS  into  the  IP,  R  represents  rotation  around  the 
IP  through  a  Unite  angle  y,  and  IT-1  takes  the  pole  of  the  TRS  into  the  IP. 

/  cosy  siny  0\ 

R  =  —siny  cosy  0  ,  (16) 

\o  01/ 

Z1 
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where  (X,Y)  denote  the  coordinates  of  the  IP  in  the  CRS,  and  x  is  the  finite 
angle  of  rotation  around  the  IP.  W  has  the  same  form  as  S  except  for  the 
replacements  X  — ►  — £,  Y  —r]  where  (£,  r /)  are  the  coordinates  of  the  IP  in  the 
TRS,  related  to  the  familiar  polar  motion  parameters  by 

£  =  xP,r)=-yP-  (ld) 

(Only  the  part  of  W  that  is  linear  in  (£,77)  is  retained  in  VLBI  analysis  at 
present,  the  higher  order  part  being  deemed  negligible.)  S  represents  precession 
and  nutation,  and  IT,  the  “polar  motion”.  X,Y,  £,77,  and  x  constitute  the  Earth 
orientation  parameters.  The  definition  of  the  CEP  envisages  that  the  spectra  of 
the  first  four  of  these  contain  only  low  frequency  components.  For  instance, 

X  =  ^[xLr  sin  0T  +  A't  cos  9t\,  Y  =  ^[1?T  cos  0T  +  B'T  sin  0T\,  (2a) 

T  T 

0T  =  rfi0f  +  ay,  (26) 

where  Oq  is  the  mean  sidereal  rotation  rate  of  the  Earth,  and  the  sum  is  over 
spectral  components  with  0  <  r  <  (1/2). 

In  current  practice,  the  expression  used  for  S  is  in  terms  of  the  nutations 
in  longitude  and  obliquity,  A  ip  and  Ac,  not  the  Cartesian  (X,  Y);  the  transfor¬ 
mation  between  these  may  be  found  in  Capitaine  (1990). 

One  sees  readily  that  the  coordinates  (Xj,Yj)  of  the  terrestial  pole  (pole 
of  the  TRS)  in  the  CRS  are  just  the  (1,3)  and  (2,3)  elements  of  the  matrix  T-1, 
and  hence  that 

Xj  ~  X  —  £  cos  x  +  77  sin  y,  YT  =  Y  -  £  sin  x  -  77  cos  y,  (3) 

to  the  first  order  in  (X,Y)  and  (£,77).  The  largest  of  the  higher  order  terms  in 
Xj  is  XF ,  and  that  in  Y  is  YF,  where 

F  =  [X(£cosy  -  ??sinx)  +  T(£sinx  +  77 cos x)]-  (4) 

The  magnitude  of  F  is  quite  small.  Precession  at  the  rate  p  ~  50"  per  year 
contributes  lOOpsineo  ~  2000"  ~  0.01  radian  to  X  over  a  century  while  £  and 
77  are  under  0"4  ~  2  X  10-6  radian.  So  F  does  not  exceed  2  X  10-8  over  a 
century,  the  corresponding  correction  to  Xj  of  (3)  would  be  less  than  40  pas 
over  a  century.  The  periodic  part  of  the  correction  is  under  0.2  pas  for  both  X 
and  Y.  The  linear  approximation  (3)  is  therefore  sufficient  for  present  purposes. 
The  expression  for  the  complex  combination  Xj  =  (Xj  +  iYj)  is  more  revealing 
than  eqs.  (3)  as  they  stand: 

XT  =  X  —  £e*'x,  (X  =  X  +  iY,  £  =  £  +  ;? 7).  (5) 

The  expansions  (2)  for  X  and  Y  lead  to 

X  =  Y^CTetTQot,  (6) 

T 

where  the  sum  is  now  over  both  positive  and  negative  (prograde  and  retrograde) 
frequencies  in  the  range  (  —  flo/2)  to  Oq/2.  The  expansion  of  £  is  similar. 
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The  nature  of  the  terms  containing  the  polar  motion  variables  in  eqs.  (3) 
or  (5)  deserves  special  attention:  with  £  and  r]  taken  (as  in  routine  VLBI  anal¬ 
ysis)  to  contain  no  high  frequency  components,  these  terms  represent  purely 
prograde  diurnal  components  of  the  motion  of  the  terrestrial  pole  in  space  if  we 
leave  aside  for  the  moment  the  extremely  small  terms  in  other  frequency  bands 
(to  be  considered  later)  that  arise  from  the  diurnal  and  semidiurnal  UT1  vari¬ 
ations  contained  in  y.  There  is  no  place  within  this  scheme  for  semidiurnal  or 
retrograde  diurnal  signals  in  X  and  Y  due,  for  instance,  to  ocean  or  atmospheric 
tides.  Herring  and  Dong  (1994)  overcame  this  limitation  by  a  scheme  which  is 
equivalent  to  replacing  the  slowly  varying  £  in  (5)  by 

£o  +  ^2  £neinQot ,  (7) 

wherein  £n  for  every  n  contains  only  low  frequencies  (between  (  —  Do/2)  and 
(Do/2).)  When  this  is  introduced  in  (5),  £o  takes  the  place  of  £,  while  the  term 
for  any  n  ^  0  represents  signals  within  a  band  of  width  Do  centered  at  the 
frequency  (n  +  1)Dq  in  space. 


3.  Is  an  Intermediate  Pole  necessary? 

Long  usage  of  the  factorization  (1)  of  the  transformation  between  the  CRS  and 
the  TRS,  and  the  accompanying  identification  of  the  polar  motion  parameters 
in  W  with  the  coordinates  of  an  Intermediate  Pole  (hitherto  the  CEP)  has  led 
to  a  universal  and  unquestioning  belief  that  an  Intermediate  Pole  is  essential  to 
make  it  possible  to  estimate  both  the  precession-nutation  and  the  geophysically 
induced  rotation  variations  (such  as  the  Chandler  wobble)  as  low  frequency 
signals  from  VLBI  data. 

Remember,  however,  that  the  VLBI  observables  depend  on  Earth  orienta¬ 
tion  only  through  three  independent  parameters,  say  ,  and  the  angle  of 

rotation  yy  around  the  terrestrial  pole.  For  the  estimation  of  £  and  r]  from  the 
data,  the  manner  in  which  they  enter  into  the  prograde  diurnal  part  of  (Xy,  ly) 
is  all  that  is  relevant,  not  any  interpretation  as  coordinates  of  some  intermediate 
pole.  Independently  of  this  line  of  reasoning,  one  of  us  (Mathews,  1999)  had 
proposed  a  scheme  for  the  modeling  and  estimation  of  signals  in  any  frequency 
band  (n  —  l/2)Do  to  (n  +  l/2)Do  of  interest  in  the  CRS  or  the  TRS.  If  the  ter¬ 
restrial  pole  itself  (rather  than  some  Intemediate  Pole)  is  chosen  as  the  reference 
pole,  then  the  scheme  reduces  to  (i)  expanding  (Xy,ly)  as 

XT(t)  =  X0(t)  +  ^pW(t)  cos (nD0/)  -  Yn(t)  sin(nD0/)], 

YT(t)  =  Y0(t)  +  J2[Xn  ( t )  sin(raDot)  +  Yn(t)  cos(nD0/)],  (8) 

with  the  spectrum  of  every  Xn(t)  and  Yn(t)  confined  to  low  frequencies  (0  to 
Do/2),  and  (ii)  estimating  as  many  pairs  (Xn(t),Yn(t))  as  are  desired,  subject 
to  limitations  imposed  by  the  quality  of  the  data.  The  terms  corresponding  to  a 
given  n  /  0  in  Xy(/)  and  Yj(t)  would  contain  frequencies  between  (n  —  l/2)Do 
and  (n  +  1/2)Dq  in  the  CRS,  prograde  (retrograde)  for  positive  (negative)  n; 
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and  (Xo(t),Yo(t))  represents  the  low  frequency  band  (  — l/20o  to  l/20o)  in  the 
CRS.  This  is  manifest  from  the  complex  version  of  (8): 

XT(t)  =  10  (/)  +  Xn(t)emU°\  (9a) 

the  spectral  expansion  of  Xn(t)  being 

Xn(t)  =  Y,  Xn^ot.  (9b) 

—  1/2<t  <1/2 

It  is  now  evident  that  if  (8)  were  truncated  by  dropping  all  but  the  n  =  1 
term  in  the  sum,  the  result  would  be  just  the  same  as  (3)  with  only  one  difference: 
the  replacement  of  y  in  (3)  by  fled  in  (8).  It  may  be  verified  that  the  contribution 
to  Xj  or  Yj  in  (3)  from  the  largest  of  the  diurnal  and  semidiurnal  terms  in 
the  difference  (y  —  fled)  is  under  10-6  mas  and  hence  entirely  negligible.  The 
low  frequency  part  of  (y  —  0(d)  would  contribute  to  the  low  frequency  part  of 
(Xt,Xt)-  The  difference  between  (Xo,Yo)  and  (X,Y)  of  eq.  (3)  that  arises 
therefrom  is  well  under  the  gas  level,  except  for  the  contribution  from  the  18.6- 
year  term  in  y  which  is  still  under  5  gas.  Except  for  this  tiny  difference,  (Xo,  To) 
of  (8)  is  just  the  low  frequency  (X,  Y)  of  (3).  More  importantly,  (Xi,Yi)  is 
seen  to  be  just  the  same  as  (£,77),  though  the  former  is  not  identified  with  the 
terrestrial  motion  of  any  Intermediate  Pole:  no  such  pole  was  employed  in  the 
derivation  of  (8). 

We  may  therefore  draw  from  the  correspondence  above  the  important  con¬ 
clusion  that  the  invocation  of  an  Intermediate  Pole  in  the  traditional  decompo¬ 
sition  of  the  transformation  T  is  not  of  any  real  relevance  to  the  definition  of 
the  Earth  orientation  parameters,  and  that  one  may  work  just  as  well  with  the 
terrestrial  pole  itself,  and  the  parameters  (Xo,Yo),  (Xi,Yi),  and  y y  associated 
with  it. 

A  further  important  fact  is  that  (8)  already  contains  additional  terms  rep¬ 
resenting  frequency  bands  other  than  the  two  present  in  (3):  n  =  2,3,- ••  for 
prograde  semidiurnal,  terdiurnal,-  •  •  bands  of  frequencies  in  space,  and  n  = 
—  1,  —  2,  •  •  •  for  retrograde  diurnal,  semidiurnal,-  •  •  ones.  The  signals  in  one  or 
more  of  these  bands  (due,  for  instance,  to  ocean  or  atmospheric  tides)  can  be 
found  simply  by  estimating  the  corresponding  low  frequency  functions  (Xn,Yn) 
along  with  those  for  n  =  0  and  1  and  y y.  This  procedure  is  essentially  equiva¬ 
lent  to,  but  more  direct  than,  that  pioneered  by  Herring  and  Dong  (1994);  it  has 
been  found  successful  in  tests  done  recently  by  Ch.  Bizouard,  A.-M.  Gontier, 
and  N.  Capitaine  (private  communication,  2000). 

Yet  another  fact  of  some  importance  is  that,  there  being  no  intermediate 
pole,  there  is  no  scope  for  the  oft-debated  question  as  to  whether  signals  in  one 
frequeny  band  or  the  other  are  to  be  visualized  as  and  estimated  in  the  CRS  or 
in  the  TRS,  or  for  the  related  question  as  to  whether  signals  arising  from  some 
particular  physical  cause  (e.g.,  the  diurnal  and  semidiurnal  nutations  produced 
by  external  gravitation)  should  be  modeled  in  the  CRS  or  the  TRS.  All  terms  in 
expressions  (8)  and  (9)  are,  unambiguously,  parts  of  the  motion  of  the  terrestrial 
pole  in  the  CRS. 
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4.  Modeling  and  estimation 

Modeling  is  done  by  using  expressions  of  the  form  (2)  for  each  of  the  Xn(t)  and 
Yn(t ),  with  values  derived  from  theory  for  the  coefficients  AntT ,  A'n  T,  BUtT,  B'n  T. 
For  nutations  caused  by  external  gravitating  bodies,  one  has,  besides  the  main 
terms  corresponding  to  n  =  0,  also  the  recently  discovered  terms  belonging  to 
n  =  1  and  2.  For  variations  of  geophysical  origin,  signals  having  terrestrial 
frequencies  <rflo  with  a  in  the  range  ( r  —  1/2)  <  a  <  (r  +  1/2)  are  represented 
by  the  term  with  n=(r+l)  in  (9).  In  particular,  one  has  a  term  Xr_|_ijTe*(,’+1)fiot 
corresponding  to  a  wobble  of  frequency  cOo  and  amplitude  fh(a ),  (a  =  r  +  r): 

fh(a)  .  . 

Xr+hT  =  -——,  ( o  =  r  +  r ).  (10) 

1  +  cr 

Modeling  of  such  a  term  can  be  done  if  a  satisfactory  theoretical  model  is  avail¬ 
able  for  the  computation  of  fh(a).  A  wobble  fh(a)  due  to  a  prograde  diurnal 
ocean  tide,  for  example,  would  contribute  a  term  X2tTe^2+T^ot  in  the  semidiur¬ 
nal  band  to  X y(t);  semidiurnal  nutation  terms  excited  by  external  gravitation 
also  appear  in  the  same  band. 


5.  Kinematic  relations 

The  components  Oq^i,  om2 ,  and  Oq(1  +  m3 )  °f  the  Earth’s  angular  velocity 
vector  fi  in  the  TRS  are  the  off-diagonal  elements  of  the  antisymmetric  matrix 
M  =  T~1T,  where  the  transformation  matrix  T  is  just  RS  in  the  absence  of  any 
Intermediate  Pole.  Evaluation  of  M  leads  to  the  kinematic  relations 


flQmi  =  L2 3  =  (X  +  XA)  siny  -  (Y  +  Y A)  cosy,  (11a) 

<fl0m2  =  L31  =  (X  +  XA)  cosy  +  (Y  +  Y A)  siny,  (116) 

fio(l  T  m3)  =  T12  =  fi  +  A3,  (Hc) 

A  =  XX  +  YY,  n  =  x,  A3  =  (1/2)(YX  -  XY).  (12) 


Eq.  (10)  is  the  frequency  domain  version  of  the  equation  flofh  =  iXe~tQot  that 
is  equivalent  to  the  linearized  form  (A  — ►  0,  y  fled)  of  the  pair  of  eqs.  (11a), 
(lib);  this  approximation  is  more  than  adequate  for  the  purpose  for  which  the 
use  of  (10)  was  envisaged. 


6.  Definition  of  UT1 

UT1  is  intended  to  be  proportional  to  the  angle  of  rotation  of  the  Earth  in  space. 
It  would  seem  eminently  reasonable  to  require  that  the  reference  pole  around 
which  the  rotation  is  to  be  defined  be  a  physically  meaningful  one,  such  as  the 
pole  of  the  TRS  or  the  instantaneous  rotation  pole  —  not  some  “Intermediate 
Pole”  (the  CEP  or  some  modified  version  of  it)  that,  strictly  speaking,  has 
no  real  physical  existence.  It  is  pertinent  to  point  out  in  this  context  that 
UT1  variations  of  geophysical  origin  (due  to  zonal  tidal  deformations)  must 
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necessarily  be  identified  with  variations  in  rotation  around  the  pole  of  the  TRS 
because  the  third  axis  of  the  TRS  is  the  reference  axis  for  the  definition  of 
the  zonal  tides.  The  choice  of  a  different  reference  pole  for  the  fundamental 
definition  of  UT1  would  then  seem  incongruous. 

When  the  Earth  orientation  parameters  are  defined  in  relation  to  an  IP,  the 
third  of  the  kinematical  relations  has  the  general  form 

fl0(l  +  m3)  =  L\2  ~  (1/2 ){ypxp  -  xpyp)  -  [(l/2)L12(x2p  +  yp)~  (L23  xp  -  L31  yp)\, 

(13) 

instead  of  (11c),  where  the  Lij  are  defined  as  in  eqs.  (11).  Recalling  that 
L12  =  dx/dt  +  A3,  we  find  the  angle  of  rotation  around  the  pole  of  the  TRS  to 
be 


/  f2o(l  +  ^3)  —  X  +  (1/2) 


XY)dt-{  1/2)  (ypxp-xpyp)dt+  -  ■  (14) 


wherein  the  dots  represent  the  integral  of  the  square-bracketed  part  of  (13),  and 
we  have  used  the  expression  for  A3  from  (12).  One  recognizes  the  terms  following 
X  in  (14)  to  be  of  just  the  same  forms  as  the  quantities  s  and  s!  involved  in  the 
implementation  of  the  nonrotating  origin  by  Capitaine  (1990);  x  with  these 
terms  added  constitutes  the  basis  for  her  definition  of  UT1,  when  the  pole  of 
reference  is  the  CEP  or  a  redefined  IP.  If  the  pole  of  reference  were  the  pole  of 
the  TRS  itself,  the  replacements  x  XT ,  (X,Y)  — ►  (Xo,lfi)  and  (xp,—yp)  ~ 
(Xi,Y\)  would  lead  to  the  corresponding  quantity,  xx  +  (l/2)  J(YoXo  —  XoYo)dt+ 
(1/2)  f(YiX  1  —  X{Yi)dt  as  the  appropriate  basis  for  the  definition  of  UT1. 

It  would  seem  therefore  that  there  is  no  need  for  an  Intermediate  Pole  for 
the  definition  of  UT1  either. 
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